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Apologia 


Dittoed  manuscripts  are  by  definition  preliminary 0  This  one 
is  especially  mp  and  I  am  under  r©  illusions  as  to  the  rigorouaneos 
of  tbs  proofs  in  Sections  4  sod  5o  On  the  other  hand,  the  topic 
dealt  with  here  has  to  do  fundamentally  uith  approximations,  and  I 
have  a  notion  that  in'  this  instance  heuristic  is  of  more  importance 
than  rigor  o  Moreover,  if  anyone  tdshes  to  make  matters  more  rigorous 
the  proofs  provided  here  do,  1  believe,  mark  a  dear  route;  and  patience 
and  a  liberal  sprinkling  of  epsilons  will  fill  in  the  details, 

Actoowledflaentt  I  have  had  numerous  very  enlightening  discussions 
with  Allen  Newell  on  the  subject  of  this  paper. 


1.  MroductioR 

The  mood  to  deal  xjith  many  problems  of  economic  theory  in  terns 
of  aggregates  1b  obvious  0  The  complete  Walrasian  system,  and  the  more 
modern  dynamic  embellishments  of  it  are  relatively  barren  of  results  for 
macroeconomics  and  economic  policy.  Hence,  by  sheer  necessity,  we  use 
highly  aggregated  systems,  often  xdthout  having  much  more  than  the  same 
necessity  as  a  justification  for  their  use.  Perhaps  the  most  important 
result  to  date  in  the  attempt  to  justify  aggregation  under  certain 
circumstances  is  the  Longe-Hiclcs  condition,  about  which  we  shall  say 
more  later. 

The  development  of  the  Leontief  input-output  model,  and  the 
concern  with  actual,  numerical  coefficients  in  rather  highly  disaggre~ 
gated  versions  of  this  model  have  renewed  interest  in  aggregation. 
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There  was  at  first,  perhaps,  the  food  hops  that  modern  oowpctlng 
squipnent  ^xjuld  handle  matrices  of  any  3tse  we  night  be  likely  to 
concern  ourselves  vdth,  and  hence  would  obviate  the  need  for  aggrega¬ 
tion  <,  By  now,  it  is  dear  that  the  sine  of  the  models  we  should  like 
to  handle  has  outstripped  the  developraent  of  ooopgfcero  -and  is  lilcely 


to  continue  to  do  so  as  long  as  tho  time  required  to  invert  a  matrix 
increases  with  the  cube  of  the  number  of  rows  and  oolxmnso  A  nunber 
of  exporlments  have  already  boen  aade  (e.g.,  Vft&tin  and  tbrgenstem)  to 
determine  whether  aggregation  can  be  used  to  obtain  approximate  inverses 
more  economically. 


Most  discussions  of  aggregation  have  taken  place  in  the  oontext 
of  ea  algebraic  model.  The  problem  nay  be  stated  In  general  tense  as 
follouet  Tib  hare  two  eats  (not  neoeesnrily  distinct)  of  variable* ■ 

a  system  of  equations  giving  the  y*s  as  functions  of  tho 
x*s.  Tie  wish  to  know  under  ximt  drmsastanoes  there  exists  a  function, 
of  the  x*s  and  another  function,  T(y^a0jyn),  of  the  y*s, 
such  that  a  relation,  7  «£(X) ,  between  7  and  X  can  be  bieduoed  from  the 
given  system  of  equations  relating  the  x*e  and  y*sa  Sometimes  additional 
conditions  are  imposed—*. g»,  the  functions  7  and  X  are  given,  or  it  is 
required  that  the  relation  jg,  poesess  certain  properties  0  (For  a  compre¬ 
hensive  treatment,  see  Theil)« 

The  conditions  that  must  be  satisfied  in  order  for  aggregation  to 


be  strictly  possible  are  vary  severe.  Whether  these  conditions  are 
exactly  satisfied  in  any  practical  situation  Is,  of  ooursa,  not  an 
important  question.  TJe  would  be  perfectly  satisfied  with  aggregative 


models  that  gave  us  only  approximate  results— we  have  no  illusions.) 
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after  all,  that  any  modal  :w  might  aaipJuoy  is  raora  than  on  approximate 
description  of  reality.  Our  objective  in  exploring  the  aggregation 
problem  should  bo  to  sook  rules  and  c riteria — ox}.ct  or  heuriotio— that 
xjill  give  us  clues  aa  to  what  variables  to  aggregate,  and  t.111  indicate 
to  us  wader  what  oircurastanceo  aggregation  is  likely  to  yield  satisfactory 
apr  roxi  nations . 

The  lango-ffiLcks  condition  is  a  criterion  of  this  kind.  It  states 
t/jut  if  two  or  more  variables  always  move  together,  then  they  my  be 
aggregated  into  a  single  variable,  which  will  be  an  appropriately  weighed 
average  of  the  original  variables.  This  is  a  useful  criterion,  since  it 
tell.  1  us  that  we  my  aggregate  classes  of  cocnodities  that  are  perfect 
substitutes— or  that  are  approximately  so. 

At  another  level,  tbs  Lange-Hleks  condition  is  unsatisfactory, 
for  it  requires  that  we  know  as  a  datum  which  variables  neve  together  0 
In  actual  fact,  it  may  be  part  of  our  problem  to  discover  this.  He  may 
be  confronted  with  a  dynamical  system  involving  a  large  number  of  varia¬ 
bles,  and  may  have  to  infer  from  the  equations  of  the  system  which  varia¬ 
bles  will  move  together— or  will  behave  nearly  enough  in  this  manner  to 
warrant  aggregation.  This  is  the  essential  problem  we  shall  set  our¬ 
selves  hero:  to  determine  conditions  that,  if  satisfied  by  the  equating 
of  a  (linear)  dynamical  system,  -.dll  pewit  aPOTOodmate  aggregation  of 
variables.  Note  that  we  will  be  interested  in  sufficient,  rather  than 
necessary,  conditions.  Hence,  we  my  also  view  our  task  as  that  of 
discovering  one  or  more  classes  of  dynamioal  systems  that  permit  aggrega¬ 
tion  of  their  variables 0 
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20  Properties  of  "Nc&rly-Decoapooable"  Matrices 


For  purposes  of  exposition,  it  will  be  oonveniont  to  lirdt  our¬ 
selves,  for  the  nest  part,  to  stochastic  natrices— that  is,  am  matrices 
all  of  whose  elements  ore  non-negative  and  with  the  om  of  the  elements 
in  each  row  equal  to  It  Pi;j>0  (l,j  -  l,..4,n),  awl  •  1  (j  •  l,...,n). 
Hit  insults  nay  readily  be  generalised  (Robert  Solow,  "On  the  Structure  of 
linear  Models,"  Econometrics,  January,  1952,  pp.  42-43) ,  but  the  restriction 
will  perhaps  make  the  intuitive  basis  of  the  analysis  dearer  than  would  a 
more  general  treatment.  VflLth  this  interpretation,  the  index  i  runs  over 
the  £  possible  states  of  the  system,  p^(t)  represents  the  probability 
that  the  system  is  in  the  state  at  time  t,  and  p  represents  the 
conditional  probability  that  the  system,  if  in  state  1  at  tine  t,  will  be 
in  state  J  at  time  (t+1). 


Mow  suppose  that  by  an  approximate  permutation  of  rows  and  oorres— 


ponding  columns  the  stochastic  matrix 
lows: 


"INI 


can  be  arranged  at  fol- 


(2.1) 


P» 


P« 

A‘. 


P* 

1 


U 


whore  the  P^  are  square  sutnatrlcos  and  the  remaining  elements,  not 
displayed,  ore  all  zero.  Then  the  matrix  is  completely  decomposable. 
Suppose  further  that  no  no  of  the  submatrices,  P« ,  are  themselves 
decomposable. 

Let  the  system,  at  time  V  be  in  the  state  X^,  where  1^  is  one 
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of  tbo  states  belonging  to  the  suboatrix  P’ ,  Then  it  is  obvious 
that  at  all  subsequent  tinoo  the  ryotem  will  continue  to  be  in  one  of 
the  states  belonging  to  P’.— that  the  probability  of  transition  to  any 
state  outside  this  subsot  is  zero*  Further,  if  the  initial  position  of 
the  system  is  given  by  a  probability  distribution  over  the  states  belonging 
to  P£,  all  other  states  having  zero  probability,  then  the  limiting  distri* 
button  of  the  system,  after  the  passage  of  an  indefinite  period  of  time, 
mill  be  independent  of  this  Initial  distribution*  If,  on  the  other  hand, 
the  initial  distribution  assigns  non-zero  probabilities  to  states  belonging 
to  more  than  one  of  the  indecomposable  subraatrices,  P^4  then  the  limiting 
distribution  will  not  be  independent  of  the  initial  distribution* 

Let  the  initial  distribution  be  p^(0),  and  define* 

(2.2)  Pi/I  -  Pi/pj 
then  the  final  distribution  uill  bet 

(2.3)  PA  "  Pl(°)*Pi/i 

where  the  (but  not  the  p^)  are  independent  of  the  initial  distribution, 
The  Pjy^  are  the  characteristic  vectors  of  the  suboatrlx  P£  corresponding 
to  the  largest  characteristic  root,  X^**  1»  Further,  the  distribution  of 
(2*3)  is  an  equilibrium  distribution  in  that,  if  this  distribution  is 
once  attained  by  the  system,  it  uill  not  change  over  time. 

Now,  suppose  the  stochastic  matrix  to  be  slightly  altered  so  that 
there  is  now  at  least  one  path  with  non-zero  probability  leading,  in 
not  more  than  £  steps,  from  any  one  state  of  the  system  to  ary  other 
state*  That  is  to  say,  ue  suppose  certain  of  the  zero  elements  in  P1 
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that  lie  outside  the  autmatrloee,  P^,  to  be  replaced  tor  very  email 
non-eero  elements,  and  certain  elements  Inside  the  P|  to  be  corres¬ 
pondingly  reduced  eo  as  to  preserve  the  neen  condition.  We  will 
refer  to  the  modified  matrix,  P,  obtained  from  P*  In  this  fashion,  as 
a  "imarly-decongposable"  matrix.  The  term  "nearly-deoampooable, 0  like 
the  tem  "very  snail  nan-aero  elements”.  Is  not,  of  oouroe,  precise*  We 
will  be  ooatenb  to  leave  tho  criterion  of  aaaHneee  undefined.  The  smaller 
these  dements  are,  the  closer  to  the  exact  solution  will  be  the  appraxL- 
nation  to  be  diagnosed  here,  but  we  shall  not  attseqpt  to  evaluate  explicitly 
the  goodness  of  the  approximations. 

Comparing  the  original  decomposable  matrix,  P*,  and  the  modified, 
nsarly-dsocapoeebls  natrlx,  P,  we  any  readily  reach  the  following  ooo- 
alneiona,  wULoh  depend  on  the  fact  that  the  characteristic  roots  of  a 
natrlx  are  continuous  functions  of  the  natrlx  elsmanfeat 

1*  Unity  will  be  the  largest  characteristic  root  of  both  P*  and 
P.  It  will,  be  a  root  of  multiplicity  II  of  P*  (since  sad  indecomposable 
submatrix  of  P*  will  have  unity  as  a  simple  root);  bub  a  simple  root  of  P. 

2.  The  natrlx  P*  will  have  (n-U)  roots  that  ore  less  than  1  In 
absolute  value)  tho  natrlx  P  will  bars  a  set  of  (n-tl)  roots  In  one-one 
correspondence  to  these— the  differenoe  between  corresponding  roots 
of  P*  and  P  being  a  continuous  function  of  the  differences  between  the 
too  coefficients  of  the  oatrloee. 

3.  The  matrix  P  will  have  (N-i)  roots  very  close  to  unity.  Us 
will  suppose  the  differenoe  between  P*  and  P  to  be  sufficiently  small 
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that  the  difference  between  unity  and  the  enadloct  of  these  roots  iu 
■very  ouch  smaller  («ug.,  of  the  order  of  one  percwnt)  than  the  difference 
between  unity  and  the  largest  of  the  refraining  roots  of  P„ 

This  isomorphism  betvreon  the  characteristic  roots  of  P*  and  the 
roots  of  P  allows  us  to  nake  a  qualitative  analysis  of  the  structure 
of  the  latter  system*,  The  roots  of  P  belong  to  tliree  classes: 

(l)  the  simple  root  at  unity;  (2)  the  (B-l)  next  largest  roots  (cor¬ 
responding  to  tho  remaining  roots  of  P*  at  unity);  and  (3)  the  (n-N) 
smaller  roots.  The  latter  roots  can,  in  turn,  be  partitioned  into  sub¬ 
sets  corresponding  to  the  subaatricea  of  P1 » 

Ue  see  that  the  roots  of  class  (3)  are  involved  in  the  process 
whereby  equilibrium  is  reached  in  the  distribution  within  each  of  the 
Indecomposable  subsets  of  states  of  the  system.  The  roots  of  class  (2) 
are  involved  in  the  process  of  attaining  equilibrium  among  subsets.  The 
largest  root  is  associated  vrith  the  steady-state  distribution  of  the 
entire  system. 

The  smaller  a  characteristic  root,  the  more  quickly  is  the 

transient  damped  with  which  this  root  is  associated.  Vfe  have  assumed 

the  difference  betvreon  Pf  and  P  to  be  so  small  that  the  roots  of  class  (3) 

will  be  smaller  than  the  roots  of  class  (2).  Bence,  equilibrium  within 
(the  equilibrium  among  subsets.  In  first  approximation,  we  may  view) 
uWQoets  will  be  attained  more  quickly  lihaE^ne  equilibrating  process  as 

proceeding  in  two  successive  steps: 

(1)  wa  reach  tha  equilibrium  represented  by  (2.3),  corresponding 
to  the  equilibrliza  of  the  completely  decomposable  system  for  the  same 
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Initial  distribution* 

(2)  maintaining  this  equilibrium  vdthin  each  subsot  «.+,  all 
times,  mb  attain,  through  a  slow  change  in  the  p^.,  the  general  aquili** 
brim  of  the  angina  system,, 

Sotlee  that  in  the  seoond  stage  of  this  process,  all  the 
within  any  single  subset  vary  proportionately— hence  satisfy  the  Hicks- 
Tenge  condition  for  aggregation*  That  1b,  the  long-range  dynamic  behavior 
of  the  system  depends  only  on  the  p^a  As  for  the  short-range  dynamic 
behavior  of  the  system,  mb  may  regard  it  as  composed  of  N  independent 
indecomposable  parts— corresponding  to  the  submat rice 3  P^— and  may 
solve  eeoh  of  these  subsystems  independently  of  the  others. 


Let  Bj.  be  the  number  of  rows  (and  column)  in  P^,  and  heme 
the  ntnber  of  states  in  the  corresponding  subset  (Enj  ■  n).  Then,  by 
the  process  of  decomposition  just  described,  ue  have  replaced  a  dynamical 

process  represented  by  a  zpm  matrix  uith  a  superposition  of  dynamical 

.»  -  - 

Frocosses  represented  by  H  matrices  of  siae  .a.,  respectively* 

and  one  additional  matrix  of  siae  |bg.  Later,  we  nha^l  show  that  tMa 

4fWBP^^^a-peyait  s,  U£.  &> ..replay  S&&  JfaKfe  of  iBOgldBS  ^ttJBfadaL  T 
uith  the  much  less  laborious  task  of  lnrertinr  the  matrices  of  the 
teWWPPPPd  grates* 

3o  A  Physical  Illustration 


Before  we  proceed  uith  a  more  careful  statement  of  the  mathe¬ 
matics  that  underlies  our  coialysia,  it  may  be  useful  to  provide  an 
example  of  a  physical  system  that  can  be,  approximately,  decomposed 
in  the  manner  just  escribed.  The  Illustration  will  be  useful,  also, 
in  showing  that  our  conclusions  do  not  depend  on  the  stochastic  proper¬ 
ties  of  the  matrix  ue  have  discussed  thus  far.  We  shall  see  that  the 
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principle  of  aggregation  t,o  are  employing  ia  essentially  that  which 
justifies  tho  replacement  of  micro  var.'  aMos  by  macro  variables  in 
classical  thonaodyiinmice  „ 


Lot  ns  consider  an  iraagi’iary  building,  wboea  outa3.de  '.rails  pro¬ 
vide  perfect  thermal  insulation  from  tho  omdrorosiifc*  The  building  is 
divided  into  a  large  number  of  rooms,  the  walls  between  them  being  good* 
but  not  perfect*  insulators „  Pa  oh  room  is  divided  into  a  number  of  offices 
by  partitions,,  The  partitions  are  only  poor  insulators*  A  thermcooter 
hangs  in  each  of  the  offices. 

Let  ud  suppose  that  at  time  te  the  various  offices  within  the 
building  are  in  a  state  of  thermal  disequilibrium — there  is  a  vrf.de 


variation  in  temperature  from  Office  to  office  and  from  room  to  room* 

If  m  take  nmr  temperature  readings  at  time  t^,  several  hours  after  t^, 
what  trill  we  find?  At  t^,  there  will  be  very  little  variation  in  tem¬ 
perature  among  tho  offices  that  are  in  oach  single  roam*  bat  thero  ray 
still  be  large  temperature  variations  between  rooms.  If  we  take  readings 
again  at  time  t^9  several  days  after  t^,  we  may  find  that  there  is  an 
almost  uniform  temperature  throughout  the  building*  tho  temperature 


differences  between  rooms  having  virtually  disappeared  u 

The  uell-Iaaoun  equations  for  the  diffusion  of  hoat  allow  us  to 
represent  this  situation  by  a  system  of  differential  equations-- or 
approximately  by  a  system  of  difference  equations. 

Let  TL.(t)  be  tho  temperature  of  the  1th  office  \Moh  is  in  the 
room*  at  time  t*  Let  T(t)  be  the  vector  consisting  of  these  tempera- 
a.  oomponanta,  T(t)  -  . ^  ^  ^  T^} 
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Then  vm  ’.All  haver 
(3.1)  T(t+1)  -  A*T(t) 

wharo  A  la  a  matrix  whose  coefficients,  a^,  represent  the  rate  a  of 
boat  transfer  between  office  i  and  office  j  per  degree  difference  in 
temperature  s 

Ue  idsh  non  to  represent  ths  fact  that  a  temperature  equilibrium 
vAthln  each  room  is  reached  rathor  rapidly,  ’Hie  a  temperature  equilibrium 
among  roomo  io  reached  only  slowly,  This  tAll  be  the  caso  if  the  a 
are  generally  large  when  i  and  j  are  offices  in  the  some  room,  and  if 
the  a^  are  close  to  zero  when  i  and  j  are  offices  in  different  rooms— 
that  is  to  cay,  if  the  matrix  A  is  nearly  decomposable. 

lien  this  is  the  case,  and  as  lox^j  as  ue  are  not  inbereoted  in 
the  high  frequency  fluctuations  in  temperature  among  offloee  in  the 
same  room,  tre  can  learn  all  tie  want  to  !a»w  about  the  dynamics  of  this 
system  by  placing  a  single  thermometer  In  each  room-— it  is  unnecessary 
to  place  a  thermometer  in  each  office. 

4.  Time  Paths  of  Ilearly  Decomposable  Systems 


Vfc  will  return  now  to  the  special  case  of  stochastic  systems 
and  examine  the  dynamical  process  in  more  detail.  The  system  is  given  byt 


n 


(4.1)  (j-l,..„n) 

If  we  solve  thic  system  of  difference  equations  to  obtain  the 


p^e  as  explicit  functions  of  time,  we  get* 

M  piw  ‘S?i(  \ 


(J“l,...,n) 
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But  we  know  also  that* 

a.3)  p4(»v 

Nov  it  ie  possible  (proridod  that  all  tho  character! otic  ,4  # 
are  distinct)  to  represent  the  right-band  rise  of  (4*3)  as  a  sun  that 
is  equal  tern  by  tom  to  the  right-hand  side  of  (4*2).  (see  tfeddttxtrn, 
pp.  25-6)  To  do  this  tie  proceed  as  follows: 


For  angr  mm  natxd xt  ^with  g  distinct  chareeterlstie  roots,  we 

oaa  find  a  set  of  a  ideapotent  na  trices,  o'*]  (f-l,..*,n),  with  the 
fallowing  propertieet 

(i)  a  «  a^;l  -  (idwapotence) 

(ii)  •  a*°*  -  0#  tor  f  fa 

(iv) 

Haw,  consider  the  natrlx  At.  Taking  the  t**1  power  of  (iv) 

and  using  (i)  and  (ii),  we  find.#,  for  A*  « ||aj^| 
t  \  (t)  a  .  fc  (f) 

(  )  f 

Using  this  representation,  we  nay  express  P  •  P*  thus  x 
/.  |\  «(*)  •  2  V*  n  (l* 

(4.4)  Pjj  j 

where,  the  tr®’1  are  aatrlee*  satisfying  (i),  (ii)  and  (iii)  above, 
substituting  (4.4)  in  (4.3)  and  ooaparing  ths  resulting  expression  with 
(4.2),  we  get! 

U*5)  *3<  ’A*01  ”i) 
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Suppose  now  that  P  a|  |p^  j  i®  nearly  do conpo sable,  and  let  P* 
be  the  corresponding  decomposable  matrix  obtained  from  P  by  an  appro¬ 
priate  limiting  process.  By  examination  of  the  limiting  process  ;re 
can  again  divide  the  roots  of  P  into  three  sets:  (l)  X^»  lj  (2) 
(corresponding  to  the  remaining  roots  of  unity  in  P*  )j  (3)  XuvLfco.opXnc 
Consider  one  of  the  indecomposable  submit  xdcos,  P^,  of  P*  and 
designate  its  roots  by  »!„ 
in  (4.4)  tie  will  have: 

(4.6)  p:t  »  ♦  *£*  yjt  no&& 

I  fatfn  * 

Next  oonsider  the  gsqj  matrices  that  are  formed  from  the  above 
matrices  by  bordering  them  uith  an  appropriate  nunbir  of  rows  and 
having  element*  mo*  ’.Js  il eel  pit »  theee  bordered 

matrices  with  tho  same  symbols  as  are  employed  for  the  corresponding 
HjX  wij  matrices  In  (4*6)  •  Zt  will  be  dear  from  tho  context  In  each 
instance  which  set  of  matrices  is  intended. 

Mow,  we  nay  write; 

(4<>7)  P*  *  P<  ♦oo.«P«  *0..«P±  -  S  n*(I)  ♦  £  X*  rt'W 

*  1  *  ' 

and  similarly  for  the  t**1  power  of  P»: 

14.8)  P**  «  X|A  tT«W*  fi  XJ*  tt«(^ 

a 

low,  for  P  e  we  will  have  tbs  corresponding  representation: 

(4o9)  P*  -  ♦  !  x£n4^  ♦  £ 

M  *  f**l  * 

We  observe  that  the  last  summation  in  (4*9)  will  be  as  nearly 
equal  as  we  please,  tem  by  tern,  to  the  corresponding  sanation  in 


o  In  exactly  the  same  uay  as 
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(4o6j  i 


provided  only  that  vo  take  P  sufficiently  close  to  P’ 


We  will  have  also  the  approximate  equality; 

UolO)  JAn»(l)^wW  ♦ 

But  in  (4ol0)  the  correspondence  cannot  bo  tara-by-term,,  far 
each  column  of  has  as  elements  the  p£P  while  each  column  of  ft  ^ 
has  as  elements  the  pi/lo  If  i  and  V  both  belong  to  the  subset  ln 
then  for  moderate  or  large  t„  we 
U«ll) 


p)<t)  Vp  , 
viy  p  ty  pi»/i 

If  i  belongs  to  I  and  J  to  J  y  le  then  we  will  hare, 


approximately 

(4ol2) 


a.  (a) 


ys . 


^  i^ci  *  *i *y*t 

and  using  (4<>11)p 

(4ol3)  "  isfeljij  VV/1  Pi’4s  PJ/J 

(t) 

*  fyjPlJ 

^  £  PiVl 

Equation  (4«12)  assumes  the  P^jj,  bo  be  very  snail  when  i1  and 
jo  p«long  to  different  subsets*  In  that  casec  a  transition  from  state 
1  to  state  J  in  t  periods  nay  be  aasuned  to  inrolve  a  single  transition 
fren  one  subset  to  another  —  paths  involving  several  such  transitions 


'f/  By  ”  moderate  or  large  t”  we  mean  t  sufficiently  large  that  X  »  0 
for  ■  MHco..,nn  y  0  for  f-lP...eN0 
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having  a  vanishingly  small  probability  ,  Hones,  on  the  right-hand 
side  of  (  4.12),  we  have  retained  only  the  terms  involving  p  t  ( 
raised  to  the  first  power.  Similarly,  for  reasonably  large  t,  the 
terms  in  which  a is  very  snail  or  very  close  to  (t-l)  will  be  negligible 
oonpored  viith  the  sum  of  terns  in  which  a  is  of  the  order  of  t/2«  Bat 
for  all  the  latter  terns,  (4.U)  will  hold,  thus  justifying  the  step 
leading  to  (4,13)* 

ften  (4*13)  vjs  eee  that  p^  depends  only  on  j,  I,  and  J,  and  Is 
independent  of  i„  But  for  noderate  or  large  t  we  also  have,  from  (4.10) 

(4.14)  Jt  ■  p  .  (  — .  » 

U  ij  Pj/J  LJ  pj 

-  *  ,  (*  -  &)'*'  Z 

'  XJ  f-2  U 


end  fron  (4.13)  end  (4,9) 


(*05)  pjj>  -  *  fy.,  -  »j  -  ^2  "if 

But  for  (4.14)  and  (4*15)  to  hold  identically  in  t,  we  nust 

hsvet 

(4.16)  tp^j  -  Pj]  (f-2,...,») 


tfc  have  now  reaclved  the  important  result  that  for  oaeh  of  the 
first  N  idonpotent  ne.trices,n  ^  ,  elensnbs  within  a  given  bloc  (e,g, 
the  bloc  of  rows  and  ooloane  V£Jj  will  vary  only  with  Pj/j*  and 
otharvdoe  be  functions  only  of  X  and  J,  This  permits  us,  for  purposes 
of  inversion,  to  ropl&oe  these  matrices  with  aggregative  matrices  having 
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a  single  row  and  column  ior  each  subset  I.  The  not  hod  of  do'ng  this 
will  be  described  In  the  next  section.. 

5.  A/yreration  and  Inversion 


Prom  (4*4),  we  see  that  if  we  can  express  a  natrix  as  a  sum 
of  idmnpotent  components,  we  can  readily  invert  it.  For,  if  we  tabs 


t  *•  -  1.,  we  obtain  ianediatelyt 

(5»1)  **  m  (i^  -  l,.«.„n) 

That  is.  to  finl  P  \  we  multiply  each  idesopotent  matrix  by  the 
reciprocal  of  the  corresponding  characteristic  root0  and  sum  over  all 
such  matrices.  Ordinarily p  this  does  not  provide  ns  with  a  useful 
aomputatibnal  procedure,,  for  the  matrices,  rr  are  difficult  to 
compute.  In  the  present  ease,  however,  we  shell  see  that  (5.1)  leads  me 
to  a  very  convenient  method  of  inversion  that  doss  not  require,  moreover, 
tbs  explicit  computation  of  the  tt 


Using  (4o7)  and  (4.9).  we  hove  inediataLyt 
(J.2)  P-1'»'[it(X)  *l2  ^  nCf)  ]  *£k  (P-)'1  i?15 

The  last  two  terms  of  (5*2)  are  found  by  inverting  the  P*. 
and  computing  the  characteristic  vectors  CtJ^J  associated  with  the 
roots  at  unity.  It  rmains  to  find  the  term  in  square  brackets. 
Suppose  we  have  any  matrix.  A  •  j  |a  _  jj  ^  such  that: 

&S) 


/  (5.3) 


Than  it  follows  that: 
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(5*4)  a(t)  -  a  ,  a(t)  (**  t) 

1J  J/J  U 

We  prove  this  by  induction: 

(5«5)  *0*  "  J*  J®|j» 

_  _  «•  _  _  B  .  ,(2) 

*i/J  J«  J'J  i/J  w 

Similarly*  if  a^j”*^  •  a^yj  we  bare: 

(5*6)  ajj*  -X#  ^  aj^  a^j  eJtJ 

"  *j/j  *1? 

la  particulart  it  follows  that  aj^^  "  *j/j  ajJ*\  ®°*  ^  A 

ba  equal  to  the  sue  in  brackets  in  (5o2)t 

(5«7)  A  -  ir(l)  ♦  S  C1 1»  ^ 

f***  r 

Than,  fron  (4.11)  „  (4.13b  and  (4.66),  wa  knew  that  A  satisfies 
(5*3) v  henoa  that: 

(5.e)  ^  -  5j/j 

where,,  as  in  (4*13) 

(5*9)  Pjj  -jEj  V,  Pj^ 

Wa  now  hare  a  practical  procedure  for  investing  the  nearly** 
decomposable  matrix,  P,* 

(X)  Wa  find  the  inverses  of  the  P*. 

(2)  Wa  compute  the  and  fora  the  aggregative  natrix  Pjj 
with  eXenents  given  by  (5«9) 

(3)  We  invert  p^r  and  find  a^1^  fron  (5.8b 

(4)  Wa  now  obtain  P  *  directly  from  (5*2) 
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The  proposed  procedure  for  matrix  inversion  can  be  illustrated 
by  a  simple  example.  Consider  the  ne arl/- =d e compos ablo  matrix: 


(6.1) 


P  - 


(6.2) 


We  take  u  F': 


pi  m 


,97^0 

*0295 

.0005 

0 

.0200 

*9900 

0 

0 

0 

0 

.9600 

.0400 

.0002 

.0002 

.0396 

.9600 

.9700 

.0300 

0 

0 

.0200 

.9800 

0 

0 

0 

0 

•9600 

•0400 

0 

0 

.0400 

.9600 

(6*3) 


(6.4)  p, 


We  readily  find: 

„  1  1*031579  -.031579 

P1  " 

-*021053  1.021053 

-1  1.043478  -.043478 

P2  - 

-oO 43478  1.043478 

We  now  ooapnte  Pjj  and  its  inverse* 

9998  .00021 

.0002  .99981 


W 


(6.5)  [jpjj  *s 


1.0002  -.0002 

•,0002  1.0002 


*1° 

lo  P. 


V 

,4  .6 

Pj/2’ 

•5  .5 
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(6C6) 


and,  tron  (5-8): 


,40008 

,.60012 

-„0Q0i 

-.0001 

,40008 

. 60012 

-.0001 

-<,0001 

,00008 

-.00012 

.5001 

,5001 

,00008 

-.00012 

,5001 

.5001 

(6.7) 


From  (6„3)  and 


— 1  A# 


(6,6),  using 

(5.2),  uo 

obtain  P  \ 

1.031659 

-.031459 

-.0001 

-.0001 

-.020973 

1.02H73 

-.0001 

-.0001 

-.00008 

-.00012 

1.043578 

-.043578 

-.00008 

-.00012 

-.043378 

1.043578 

The  sane  miaerloal  exanplo  pro  video  uo  idth  an  insight  into  the 
dynardovl  prooees  as  wanned  in  (4.9).  From  (6,1) ,  w  oayde  ^ 


(128) 2 
and  P  t 

.390089 

.579037 

.016631 

.014244 

(6.8)  P128  - 

.392503 

,586246 

.011831 

.009419 

,009465 

.01J138 

•487509 

.489888 

.011385 

.015999 

.485107 

.487509 

.  (128)2 

.200776 

.298656 

.250286 

.250282 

(6.9)  P 

.200782 

.298664 

.250279 

.250275 

.200222 

.297829 

.250973 

.250976 

0 

.200225 

.297833 

.250970 

.250973 

The  ohoracteriat 

ie  root*  of  P  are  X-l,  .9996, 

.95.  .92 

Ham  x12®  -1,  .95.  .00U.  0*  .1  x*12*5  -  X,  .001*.  0»,  0*.  That 


ie  to  sear 
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7.  Sor^s  Concluding:  Conmonta 


In  the  preceding  sections,  lie  have  analysed  the  c tincture  of 


stable  dynantc&l  systems  represented  by  nearly-daco.'apo cable  matrices  L 
Ue  have  aeon  that  such  systems  my  bo  vievod  as  composite  eyoterac, 
constructed  by  the  superposition  of:  (1)  a  number  of  highly  damped 
terms,  describing  the  dynamics  of  separate  subsystems  vrithin  the  total 
system;  and  (2)  the  remaining  terns,  lees  rapidly  damped,  describing 
the  dynamics  of  the  interconnections  among  subsystems. 

If  such  a  system  starts  from  a  position  of  di ocquilibriui,  the 
process  of  reaohing  equilibrium  may  be  divided  into  two  phases, 
corresponding  to  the  decomposition  of  the  dynamical  system  that  has 
just  been  ixslioated.  In  the  first  phase,  eaoh  of  the  safcmystems  reaches 
(or  very  dose  to)  a  short-run  equilibria:;  in  the  seoond  phase, 
the  entire  system  sms  toward  equilibrium,  each  of  the  subsystems 
retaining  vory  dose  to  its  short-run  equilibrium  throughout  the  jstooess, 
Henoo,  in  the  seoond  phase,  the  set  of  variables  in  each  mibsystms 
satisfies  (approximately)  the  Lango-Hicko  condition  and  can  therefore 


be  aggregated  into  a  single  variable. 

Bence,  the  systma  variables  in  the  oase  Just  described  can  be 
represents  as  a  two  level  hierarchy,  with  the  aggregative  variables 
at  the  higher  level,  ttou,  there  is  no  reason  'toy  ue  need  to  restrict 
ourselves  to  a  two-level  hierarchy.  For  in  such  a  hierarchy,  of 
the  subsystem  variables  at  the  lower  level  alight  be  an  aggregate  of 
variables  at  a  still  lower  level  of  aggregation.  The  natri::  of  a 
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tor—  .liwl  hierarchy,  for  ojcemple,  might  look  something  Ilk®  this: 


K  W’i| 


Q2Po 

<C  fv/ 


E ,  S. 
i  41 


.1^,  rt-  i  rv,  o.^i  | 

i.  ....  j,  ■ 

ip  O  V  j 

|-3~Aj  ,  Mj 

In  tJiis  matrix,  t  he  elements  of  the  suhaatrlces  designated  as 
ll*s  would  be  of  the  first  order  o'-‘.  smallness,  and  the  elements  of  the 
E's  and  3’e  of  toe  second  order  of  smallness.  At  too  first  level  of 
aggregation,  there  ldll  bo  four  aggregative  variables  corresponding  to 
the  four  sruhmat rices  along  too  diagonal,  respectively.  At  the  seoond 
1mA  of  aggregation,  there  Hill  bo  too  aggregative  variable®,  oorree- 
ponding  to  the  blocs  indicated  by  broken  lines. 

To  invert  a  matrix  1'lce  tho  one  depicted  above,  we  would  first 

iovort  the  matrices  then  the  two  aggregative  matrices  derived  from 

L  !  II  II 

*1j.|  ‘3  -3)  | 

rxd  ,  i  ■  i-osneetivel} , 

h  p*j  I9* 

and  finally  the  second  level,  eggreg ative  matrix. 

In  ordinary  methods  of  matrix  inversion,  the  number  of  mUlpiicp- 
tions  increases  as  the  cube  of  the  alas  of  ttu  - .  -.rax.  In  the  most 
favorable  case  of  a  two-level  hierarchy,  using  the  methods  described 

hare,  ue  will  need  to  invert  matrices  of  about  sis®  n£,  and  there  will 

* 

be  about  tf*  such  matrices  3  0  Hence,  the  mmber  of  multiplications  trill 
go  up  approximately  as  the  square  of  the  sis®  of  the  matrix.  Out  if  m, 
the  s iso  of  the  largest  matrix  to  ee  inverted  at  any  level  of  aggregation, 
roRoina  constant  as  n  increases,,  thou  the  number  of  matrices  to  be  inverted 
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will  increase  proportionately  with  n*  their  size  will  not  increase, 
and  the  total  amber  of  multiplications  ;dll  increase  only  propor¬ 
tionately  uith  n. 

It  nay  be  objected  that  decomposable  matrices  are  rare  objects* 
mathematically  speaking*  and  nsarly-doconpo sable  matrices  almost  as  much 
so.  For  if  the  elements  of  a  matrix  are  s  elected  in  any  ordinary  uay  by 
a  random  process*  the  probability  that  the  matrix  uill  be  decomposable 
is  sere.  Out  there  is  every  reason  to  believe  that  neaiMleooapomblllty 
is  a  very  ooaoon  characteristic  of  dynamical  systems  that  exist  in  the 
real  world.  Instantaneous*  or  nearly  instantaneous*  action  at  a.  distance 
is  rare*  and  sharp  boundaries  around  subsystems  relatively  omnon.  tie 
have  less  assurance  that  this  holds  for  so  dal  phenomena  than  that  it 
holds  fbr  asst  natural  phenomena*  hut  the  ftmqpmit  ooemrraaee  of  "Marly* 
diagonal*  ratrinse  In  input-output  analysis  augurs  etll  fbr  this  general 
approach  to  aggregation. 

The  notion  that  most  dynamical  systems  that  are  snoountered 
empirically  consist  of  hierarchic^  of  subsystems  —  the  linkage  within 
the  subsysteas  at  each  level  being  stronger  then  the  linkage  between  sub¬ 
systems  —  is  attractive  from  another  standpoint.  Tbs  probability  that 

a  matrix  with  dements  selected  independently  and  at  random  from  a  rec- 

* 

tangular  distribution  uill  have  stable  roots  becomes  vanishingly  small 
as  the  rise  of  the  matrix  increoeas  (Ashby).  If*  however*  the  matrix 
elements  are  selected  in  such,  a  way  that  most  matrices  ore  "nearly 
diagonal"*  stability  beoomos  an  explicable  phenomenon,  It  is  probably 
no  accident  that  most  complex  ntvvotures  ve  encounter  are  hierarohial— • 
or  at  least  have  levels. 

I  sliould  like  to  coll  attrition  to  a  number  of  discussions  in 
the  economic  literature  up>n  uhi-h  the  notion  of  nearly deowposublo 
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references  la  far  from  complete.,  since  I  did  not  originally  approval! 
this  problem  from  the  standpoint  of  aggregation,  and  became  m/are  of 
the  relation  bettreen  the  two  problems  only  at  a  relatively  late  stage 
In  my  analysis. 

1.  I  have  already  pointed  out  that  the  argument  here  may  be 
regarded  as  a  statement  of  tbe  clrovnrt^noes  under  which  the  Lango-Hiclcs 
condition  will  be  satisfied .  It  ern  an  oily  be  seoa  from  our  analysis 
that  if  the  nicrc**sycte id  *  •.>  c3yrr;;.v.i:*13y  stable,  thlr,  .rill  also  be  true 
of  the  aggregated  eyste~iT  since  *.»Uc  characteristic  rooto  of  the  aggrega¬ 
tive  matrix  are  also  roots  of  tku  criminal  matrix.  I  hove  found  that 
this  stability  theorem  had  been  proved  earlier  in  Taicntsu  Yokoyama, 

"A  Theory  of  Ooaposlte  Oeoaodi'.  ~v‘’  r.si?*.  ;  H-w-rs.  May,  1952, 

^Bkegiuk  1SSHHL  the  I4uag>»!!tifla  oowflTilen  i m  ictbUidi  and 
derives  the  stability  theorem  from  this  asmnptlom. 

2.  Z  have  never  boon  satisfied  that  Qooduin*a  Justification 
of  partial  dynamics  when  two  subsystems  are  ucdlatorally  peopled  holds 
up  if  the  coupling  is  only  fh»arlyMixrLlateml. "  On  the  basis  of  the 
present  analysis,  Z  mold  suggest  that  the  Justification  for  partial 
equilibrium  analysis  or  partial  dynamics  should  rest  on  the  property 
of  rnear-d90omposabllity , " 

3.  The  experiments  by  ’.brgenntem  and  ‘.Sdtin  with  aggregation 
in  estimating  elements  of  the  inverses  of  Leontlef  natrioee  (reported 
In  Tm>nt«0qtput  Analysis  t  an  AggnJ  vel)  appeer  t*>  have  yielded  results 

ooasiatent  \iith  the  tbsorr  c'.  .  «:  \  m-.re. 

4.  Sacuolso-i  points  o  i  tbrA  v^sregatioa  of  ooaaodities  can 
seek  its  Justification  in  oitv  v  of  two  Idndo  of  principles  that,  at 
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first  blush,  appear  rather  antithetical  to  oaoh  o uwr,  Cs  the  cac 
hand  (the  Lange-fiicks  oondltlon),  we  can  aggregate  the  parte  of  a 
subsystem  when  these  are  much  more  closely  Unhurt  with  each  other  than 
they  are  uith  the  rest  of  the  syetma.  On  the  other  hand,  \m  can 
aggregate  a  set  of  variables  if  each  of  them  is  linked  Tdth  the  remainder 
of  the  system  in  just  the  seas  uey  as  ore  the  others*  Our  analysis  of 
near^deooeposabdlity  shous  that  the  fomor  condition  is  really  a  special 
ease  of  the  latter.  For  if  1  and  i  are  rarlnbles  belonging  to  different 

eubeeta  of  a  n  wurly-deoosnposable  system,  then  p^j  is  very  small,  bet 
(t) 

piJ  »  fbr  sufficiently  large  t  is  almost  independent  of  1*  That  is  to 
say,  the  linkage  between  1  and  J  is  negligible  in  the  short  run,  and 

satisfies  the  second  condition  in  the  middle  run  for  uhioh  it  is  not 
rntjUgUOe* 

5*  Several  persons  to  uhom  Z  have  reported  orally  this  method 
of  matrix  inversion  have  oonjoctured  that  It  may  have  sons  relation  to 
the  method  proposed  by  Gabriel  Kron  in  a  number  of  his  repent  ptHlira- 
tlons.  Since  neither  I  nor  the  persons  viho  have  mentioned  this  possibility 
oan  understand  Kron's  papers,  Z  am  unable  to  ascertain  whether  this  is 
the  eaoe.  In  any  event,  he  does  not  appear  to  aaloo  explicit  use  of  the 
ohazuoteristio  roots  or  the  idempotent  oomponenta  upon  which  the  astfaod 
developed  here  depends,  nor  does  be  give  any  indication  of  tbe  oomUtAoee 
under  which  aggregation  may  be  expected  to  give  a  suf  ’latently  sleep 
cpprojdmation  to  the  inverse. 


